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Abstract 

We present an ab initio theory of core- and valence resonant inelastic x-ray scattering (RIXS) 
based on a real-space multiple scattering Green's function formalism and a quasi-boson model 
Hamiltonian. Simplifying assumptions are made which lead to an approximation of the RIXS 
spectrum in terms of a convolution of an effective x-ray absorption signal with the x-ray emission 
signal. Additional many body corrections are incorporated in terms of an effective energy dependent 
spectral function. Example calculations of RIXS are found to give qualitative agreement with 
experimental data. Our approach also yields simulations of lifetime-broadening suppressed XAS, 
as observed in high energy resolution fluorescence detection experiment (HERFD). Finally possible 
improvements to our approach are briefly discussed. 
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I. INTRODUCTION 



Resonant inelastic x-ray scattering (RIXS) is a powerful tool for probing occupied and 
unoccupied densities of states at high resolution. Moreover, the RIXS signal contains valu- 
able information about the many-body excitations of a system, e.g., those beyond the pri- 
mary quasi-particle excitation.-"- However, because the RIXS signal is described by the 
Kramers-Heisenberg equation rather than Fermi's golden rule and is sensitive to many-body 
excitations, theoretical calculations of RIXS are more difficult than those of related core- 
level spectroscopies such as x-ray absorption (XAS), x-ray emission (XES), and electron 
energy loss (EELS). Even so, models of the RIXS spectrum based on single particle band 
structure arguments can be quite useful for systems with weak electron correlations, i.e., 
sp-electron systems.-"-^ In order to account for particle-hole interactions, methods based 
on the Bethe-Salpeter equation have been employed.— In addition, there have been works 
modeling d-electron systems within the single particle picture.—"— Also, for systems with 
strong electron correlations, the Anderson impurity model and atomic multiplet theories 



have been used to explain K 



XS spectra of localized d— and /-state systems.— For more 



detailed reviews see e.g., Ref. 

In this paper we introduce a theoretical treatment of RIXS based on an approximation to 
the Kramers-Heisenberg equation which uses a real space multiple-scattering Green's func- 
tion (RSGF) formalism to describe the single particle spectrum and a quasi-boson model 
Hamiltonian to account for multi-electron (e.g. shake-up and shake-off) excitations.—*^ 
Although extensions are possible, our approach is currently limited to systems with weak 
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which also uses 



correlations. Our derivation is similar in some respects to that of Ref. 
a multiple-scattering formalism and a related treatment of inelastic losses. The main differ- 
ences are: i) our expression does not rely on a single site approximation for the single particle 
Green's function; ii) we include quasi-particle self-energy corrections based on a many-pole 
model of the dielectric function; and iii) we approximate the many-body losses via a con- 
volution with an effective spectral function that includes intrinsic and extrinsic losses and 
interference effects.— In addition, with several simplifying assumptions, we demonstrate 
that the RIXS cross section can be approximated as a convolution of XAS and XES signals. 
This simplified formula is analogous to an expression in terms of appropriate joint densi- 
ties of states.— However, our result also explicitly includes the energy dependence of the 
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dipole matrix elements. In order to calculate our approximation to the RIXS spectrum effi- 
ciently, we have implemented the theory in an extension of the real space multiple-scattering 
Green's function code FEFF9.— This RSGF technique has already been used to calculate 
several other core-level spectroscopies, including XAS, XES and EELS, as well as VIS-UV 
spectra.—"— The RSGF method has been particularly beneficial for the core-level spectro- 
scopies of complex systems, since it does not rely on periodic symmetry, which is generally 
broken by the presence of the core hole. Moreover, the approach is applicable over a broad 
range of energies. Illustrative examples are presented which yield reasonable agreement 
with available experimental RIXS. Our theory also yields simulations for related spectra, 
e.g., lifetime broadening suppressed x-ray absorption spectra, as observed in high energy 
resolution fluorescence detection (HERFD) experiments. 

The remainder of this article is ordered as follows. We begin by introducing the theory 
of RIXS based on the Kramers-Heisenberg equation and then summarize our key results. 
In particular we derive an approximate formula for the RIXS cross section in terms of a 
convolution of XES and effective XAS signals. We then present a more detailed derivation of 
RIXS in terms of quasi-particle Green's functions within the RSGF formalism. Subsequently, 
we present several illustrative calculations and compare with experimental data in a number 
of weakly correlated systems. Although our treatment is in principle more general, we 
restrict our attention in this work to systems for which the quasi-particle approximation 
is reasonable. Finally, we make a number of concluding remarks. Technical details are 
relegated to the Appendices. 

II. THEORY 

A. RIXS in terms of XAS and XES 

Below we briefly outline the basic theory of RIXS and describe the key expressions used 
in our calculations. All quantities are expressed in Hartree atomic units [e = h = m = 1) 
unless otherwise noted. Formally the resonant inelastic x-ray scattering double differential 
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cross section is given by the Kramers-Heisenberg formula^ 



(fa u \^ 



dflduj fl ^ 



Em — fl — Eq + iTM 
X 6{n-uj + Eo- Ef). (1) 



Here Q and u are the energies of the incoming and outgoing photons; Ai and A2 are the 
many-body transition operators; and {"^o), \M) and \F) are many-body electronic ground, 
intermediate, and final states with corresponding energies Eq, Em, and Ep. This formula 
for the cross section can be expressed in terms of effective one particle Green's functions 



[cf. Ref. I12I and our derivation in Sec. Oil CD ] corresponding to the intermediate and final 



many-body states 



"^'^ l^\{b\d2Q\c)\'lm \{b\d\Pg\n + E') 



dVtduj 7r Vt 

X g\n -00 + E,)g\n + E,yPdi\b)] . (2) 



Here the one-particle Green's functions g and g^ are given by 

u -h^p-Vpy + iTb 
1 



where b and c denote deep and shallow core holes, |$q) is the ground state of the valence 
electrons in the presence of core hole b, the projection operators P and Q project onto 
unoccupied or occupied states of the single particle ground state Hamiltonian, and Ef,, Ec 
are the core level energies. As shown in Appendix |X1 we can rewrite this expression in terms 
of a non-local transition operator T 

d^a ^ loo \{b\d2Q\c)\^ 
dVLdoj -K VL\uj + Eb — Ec + iF^P 

X Im [{b\T\n)g%VL -u + E^)T{fl)\b)] , (4) 

where T{Q) = [l + ^Vg'^\VL + Eb)\ Pdi. We now relate this result to the x-ray emission 
/Xe(co') and an x-ray absorption like signal /2(f2, VL — u). 

d'^cr _ 0J_ f /ie(cJi)/i(r^, Q - OJ - Ui + Eb) 
dVtdu Vt J ^ la; — cji + zFftP 



where the effective absorption coefficient /i is 

fi{n, n-u) = -hm [{b\T\n)g%n -u + E,)T{n)\b)] . (6) 

The quantity /i differs from normal x-ray absorption coefficient since the dipole transition 
operator in XAS is replaced by T{Q). If the matrix elements of g'AV are much smaller 
than unity, which is the case for all but localized excitations, we may take the leading order 
approximation, and thus relate the RIXS to the usual x-ray absorption coefficient /i(a;), i.e., 

d'^(T f fle{u}l)fJ^{^ - U - UJi + Eb) 

dVtduj VL J Icj — Wi — iFfop 

Thus we obtain a relatively simple expression for the RIXS cross section in terms of the x- 
ray absorption, x-ray emission, and a resonant denominator. Moreover, the terms in either 
expression [Eq. ([5]) or ([7])] can be calculated within the RSGF framework, as in the FEFF 
codes.— '2ii22ii2^i22 It should be noted that the above expressions (Eq. [5] and [7]) are similar 
to those given in the pioneering work of Tulkki and Aberg,— in which a derivation of 
electronic resonant Raman spectra is given in terms of multichannel scattering states, and 
applied to the K- alpha RIXS of KMn04. 



B. Multiple Scattering Theory 

We now turn our attention to the application of the multiple-scattering RSGF formalism 
to Eq. ([5]). Within this formalism the single particle Green's function can be expanded 
about the absorbing atom (see Ref. l20|) 



G{r,r\E) = - 2k 



Y,\Rl{E))Glol'o{E){Rv{E)\ 



LL' 



+ 5l,l'\Hl{E)){Rl{E)\ 



(8) 



Calculations of the RIXS cross section require both the single particle XES signal, which 
is relatively simple to calculate with FEFF9, and the effective absorption cross-section Ji. 
In order to calculate the latter we begin by rewriting Eq. ([6]) in terms of the one-electron 
density matrix p'^. 



Ji{yL,VL-u) oc {h\T\^)p''{VL~u + E^)T{Sl)\h). 



(9) 
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Note again that \b) and |c) signify states calculated in the presence of the deep or shallow 
core hole respectively. Using the fact that p = — (I/tt) Im [(7] and inserting our expression 
for the Green's function, we obtain 

LU 

pS^0L'(^-^ + ^c)](i?2'imi&), (10) 

where we have assumed that the argument Q—u+Ec is real, and used the result Im[|_R) {H\] = 
\R){R\ for energies on the real axis. Note that the energy arguments of the bras and kets 
have been omited above for the sake of brevity. We now turn to the matrix elements of the 
transition operator 

Ti,(fi) = {Rl\T{n)\b) = {Rl\ [AVg\ny + l] d,\b), (11) 

where for simplicity, we have neglected the projection operator P, i.e. approximated P=l. 
Then rewriting the Green's function g'' in spectral representation, and again inserting the 
MS expression for the Green's function in Eq. ([8]) gives 

/Oh-, 
dooi 

Y^iRlimRl) + p'llA^ - 00,)] {Rl\d,\b). (12) 

Li 

Thus in addition to the usual dipole matrix elements, we have a second term which depends 
on both energies in the problem, the incoming photon frequency Q and the energy loss Q — u. 

C. Many-Body Effects and the Quasi-Boson Model 

In this subsection we discuss the application of the quasi-boson model^ to calculations 
of inelastic loss effects in RIXS. Assuming that the absorption occurs from a deep core level 
\b) and employing the dipole approximation for the transition operators in Eq. ([1]) gives 

Ai = ^(A;|rfi|6)46 + h.c. 

k 

A2 = J2(b\d2\k)bkk + h.c.. (13) 

k 
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If we neglect exchange terms between the particle and hole, or at least assume that they 
are dealt with via an effective single particle potential, we can write the many-body ground 
state as 

\'^o) = \%)\b)\k2), (14) 

where ^2 is associated with a specific term in the sum over states in A2, \b) is the deep 
core state excited by the absorption event, and |$o) is an — 2 electron wave-function. 
Note that this approximation is only justified if k2 denotes a core electron or a high energy 
photo-electron, although we will use the approximation for valence electrons as well. This 



gives 



(15) 



where M*^'' = {k\di\b) and 

i^l^o) = ^ol^o) = {eb + ek, + K)|^o). 

Note that Ep is now the Fermi energy. Then the RIXS cross section becomes 

(Pa 1 



-Im 



unocc occ 



J2 X1^2^'(^2*^^f')*^^' 
.k\k2 k^kA 



dfldu TT Q 

(A;i|(A;3|($o|i^fc3fc4(6,6) |$o) 1^4) lA^s) 

where ^1 = i7 + Eq, ^2 = ^ + Eq — oj, and K is given by 

Kkk'{uj,uj') = G{u)clbG{uj'pCk'G{u) 



(16) 



(17) 



(18) 



and G{E) = 1/{E — H + i6) is the many-body Green's function. We now introduce a 
quasi-boson approximation to the Hamiltonian following the treatment of Ref. isl. In this 
approach the excitations of the many-body valence electronic state are represented as bosons 
while the photo-electron and hole are treated via an effective single-particle theory 



H = Hq ^ + hp + hh + Vhv + Vpv + Vph, 



(19) 
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where hh and hp are the one-particle Hamiltonians for the hole and particle respectively, 
Vhv/Vpv describes the interaction of the hole/particle with the valence electrons, 

= ^kclck; = - XI efcCfecJ., (20) 

k k 
Vpv = J2 + KkJ*<^n\ 4,C,„ (21) 

Vhv = E + KJ*^n\ CkA., (22) 

n,kik2 

and Vph describes the interaction between the photo-electron and hole. This last interaction 
term should in principle be treated via the Bethe-Salpeter equation; however, here we will 
approximate it using either a self-consistent final state rule approximation for deep core 
holes (i.e., with the screened core- hole potential of the deep core-hole), or by neglecting 
it altogether, as in the initial state rule (independent particle approximation) for valence 
holes. Experience with such models in the FEFF code shows that these approximations are 
reasonable. 

We now define |$o) ^ the ground state of the N — 2 electron system in the presence of 
the deep core hole and |$q) as the ground state of the N — 2 electron system in the 
presence of the second core hole |c) so that 

H^l^:,) = E^,\^:,); H^ = Ht^ + V;i,i, (23) 

with core level energies 

Eh — ^b — Eq^- Eq, 

E, = ec- El + E'^. (24) 

The transition matrix elements corresponding to emission (^2) may be pulled outside the 
imaginary part, and serve as an amplitude factor, i.e., 

cPa 1 oj 



dflduj TT fl 



X Im 



c 

unocc 



J2{b\4PFiE,,E^)Pd,\b) 



kik2 



(25) 



8 



Here Ei = Q + Eb, E2 = Q + Ec — uj, P is a projector onto unoccupied states of the 
ground state Hamiltonian, Q is a projector onto occupied states of the intermediate state 
Hamiltonian, 

F{E,, E,) = G\E,)G%E2)G'\E,)- (26) 

and finally, the Green's functions are calculated in the presence of the deep (6) or shallow 
(c) core hole, i.e., 

1 



G\uj) 
G^uj) 



u-iHh- El) -hi- + iFfe 



uj-{E,-E^)-hl-V^,^iY' 



(27) 



^0/ '"p " po 

Next we derive an expression for the effects of multi-electron excitations in terms of an ef- 
fective spectral function. Within the quasi-boson approximation, the following relationships 
between the eigenstates of i^o; -f^o; -^0 holci^'i^ 



|$o) = e"^''|$^); Si 



1*0) 



Aa 



E 



2 



(28) 



Here Al^*^ = V^" — VJ,^ is the difference between the intermediate and final state core hole 
potentials. If we assume only single boson excitations, we can also write 

Ar^ 



I nl 



(29) 



which will give us the correct expression to second order in the couplings when used in our 
formula for the RIXS signal. Ignoring the off-diagonal terms in V^^ we obtain 



dfldu 



1 u 
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(30) 
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Note that in the case of valence emission (valence hole), we are assuming that the core 
hole potential is negligible, hence Ec = ec, and Fc = 0. Expanding to second order in the 
amplitudes to create and annihilate bosons, and neglecting off resonant terms, gives the 
total cross section in terms of a convolution with an effective spectral function A^s. 



dVtduj 



' d^a ' 




dQdu 


sp 





(31) 



where [d'^a/dfldu]^^ is the single particle cross section as given in Eq. (E]), and the spectral 
function is given by 



A^s{Ei,E2,u;i,U2) = e ^M^M 

+ X] Wcn{E2)acn{E2)5{uJi)5{uj2 - UJn) 
n 

+ \(3tyrXEl)\'^5{uJi-UJn)5{uj2-UJn)] 



(32) 



Here Pn are the amplitudes to create or annihilate a bosonic excitation, respectively, and 
include extrinsic as well as intrinsic amplitudes (see Appendix |B]). It should be noted that 
our current formalism for the spectral function is not suited for highly correlated materials, 
although an extension of the quasi-boson model Hamiltonian is possible, as suggested in 



Ref. 



31 



and 



32 



18 and 



19 



where 



The application of the spectral function to RIXS is similar to that of Ref. 
a convolution was applied to XAS. In this paper, however, we will restrict our calculations 
to the quasiparticle approximation, i.e., with the spectral function replaced by a 5-function, 



Aes{Ei,E2,UJi,UJ2) = 6{Ui)6{u2 - 



(33) 



The use of this approximation is expected to cause the calculated spectral line-shapes to be 
more symmetric than experimental results, since the main quasiparticle peak is modeled in 
the above as a Lorentzian, while in general multi-electron excitations lead to asymmetric 
peaks so that Eq. (132|) has a Fano type main lineshape. Thus satellite peaks due to multi- 
electron excitations are also neglected. 
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FIG. 1: (color online) Calculated (left) Mn Kq, RIXS of MnO based on Eq. ([311) compared to 
experiment (right), (right). 

III. EXPERIMENT 

The experiments described here were performed at beamhne ID26 of the European Syn- 
chrotron Radiation Facihty (ESRF). The incident energy was selected by means of a pair 
of cryogenically cooled Si crystals in (311) reflection with an energy bandwidth of 0.2 eV 
(0.3 eV) at 4.9 keV (6.5 keV). The incident flux on the sample was 1 x 10^'^ photons/second 
using the fundamental peak of the undulator radiation. The beam size on the sample was 0.2 
mm vertical by 1.0 mm horizontal. Higher harmonics were suppressed by three Si mirrors 
operating in total reflection. The resonantly scattered x-rays were analyzed using the (331) 
and (400) reflection of spherically bent Ge single crystal wafers for the Ti Kp and Ka emis- 
sion, respectively. The Ge (333) reflection was used for Mn Ka. Sample, analyzer crystals 
and an avalanche photo diode were arranged in a vertical Rowland geometry {R = 1 m) at 
90 ± 3 deg scattering angle. The combined instrumental energy bandwidth was 0.8 — 1.0 
eV. All samples were purchased from Aldrich and used as is. Self-absorption effects distort 
the spectral shape and let the K absorption pre-edge region appear stronger relative to the 
edge jump. These effects are negligible in the K absorption pre-edge region and the samples 
were not diluted for the measurements. 
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FIG. 2: (color online) Calculated RIXS (left) of Ti02 based on Eq. (j3ip compared to experiment 
(right), for (top to bottom) Ti Kq,, K^, and i^vaience RIXS. 



IV. RESULTS AND DISCUSSION 



A. RIXS 



Calculations of RIXS were carried out for several materials based on the present theoret- 
ical approach using an extension of the RSGF FEFF9 code applied to Eq. fl3T|) . All results 
were calculated using self-consistent potentials and a full multiple scattering (FMS) treat- 
ment of the Green's functions for suitably large clusters centered at the core absorption site. 
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The core- hole screening was calculated using the random phase approximation (RPA),— and 
non-spherical parts of the core-hole potential were neglected. In order to simplify the calcu- 
lations the Green's functions are restricted to include only the site and angular momentum 
diagonal elements. We find that for the cases presented here, the angular momentum diag- 
onal approximation is reasonable for all but the lowest energy peaks since the overlap with 
AV is then small, and our approximation then becomes equivalent to Eq. ((Tj) in terms of XES 
and XAS where the angular momentum diagonal elements dominate due to dipole selection 
rules. To obtain better agreement with the the experimental threshold energy, we allowed a 
small shift of the calculated Fermi energy, which is typically too high by about 1 eV in the 
self-consistent FEFF9.0 calculation. In addition, overall energy shifts were added in each 
axis in order to align the calculation with the energy scale in the experiment. In the case 
of Kq, RIXS, atomic values were used for the splitting between the P1/2 and P3/2 emission 
energies.— In addition the amplitudes were taken (from simple counting arguments) to have 
a ratio ^3/2/^1/2 = 2, although this ratio is not generally accurate, since the particle-hole 
interaction mixes the hole states. Fig. [1] presents a comparison of our calculated (right) 
Mn Ka RIXS of MnO and experimental data (left). The overall agreement is qualitatively 
satisfactory: all main features of the experiment are reproduced including both the dipole 
as well as quadrupole pre-edge peaks. The main edge is also at about the correct energy, 
although the asymmetry caused by multi-electron excitations is absent in our calculation 
which is restricted to the quasiparticle level where the spectral function is given by Eq. (|33ll . 
We expect that going beyond this quasi-particle approximation for the spectral function as 
in Eq. [32] would capture some of the asymmetry since the Lorentzian spectral shape cur- 
rently used for the quasiparticle peak in these calculations would be replaced by a Fano 
type lineshape. In addition, new features could arise due to satellite peaks in the spectral 
function. The main diagonal structure in our calculation appears to be sharper that that of 
the experiment; this is possibly due to self-absorption effects in the experimental data. Note 
that the pre-edge peak for this case is basically on the diagonal; i.e., the emission energy is 
the same for the pre-edge peak as for the main edge. 

Core-hole effects are important for a variety of excited state spectoscopies, including EELS 
and XAS as well as RIXS.-*^ RIXS spectra in particular however, can give us insight into 
these effects since the intermediate and final states have different core-holes.^ In order to 
illustrate the effect of different final state core- holes, we calculated the RIXS for Ti K^, Kp, 
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and Kvaience RIXS of Ti02 Anatase. In order to obtain reasonable results for the quadrupole 
peak in the K-edge absorption, we increased the strength of the core-hole potential by using 
95% screened core-hole and 5% bare core-hole, and kept the same ratio for all core-hole 
calculations. Note that our calculation again reproduces all peaks, although the intensity 
of the calculated quadrupole peaks is weak compared to that observed in the experiment. 
There is also a noticeable effect on the spectrum due to changes in final state core-hole. 
For the RIXS, the intermediate (Is) and final (2p) core-holes are both quite localized 
and the difference AV is small, thus peaks should occur roughly on the diagonal, as seen 
by Eq. f[T^ . For the K — /3 spectrum the final state has a 3p core hole, and the core-hole 
potential has a vastly different shape than the Is core-hole potential. This causes AV to 
be large, and we expect off diagonal peaks to be present. This is indeed the case, although 
only the quadrupole peaks are off diagonal. This is due to the fact that the dipole pre- 
edge peaks are caused hj p — d hybridization between the absorbing atom and neighboring 
Ti atoms, and thus are relatively unaffected by the core-hole potential. The quadrupole 
peaks, however, are due to a direct transition to the Ti ci— states which are localized around 
the absorbing atom and are effected by the core-hole potential to a greater extent than 
the hybridized p-states. The effect is also present in the valence spectrum. In addition, 
the valence spectrum has multiple peaks due to the fact that the emission is from a broad 
valence band which is split due to solid state effects. The qualitative structure of the valence 
band is also correct in the calculation, which reproduces the double peak structure with the 
correct splitting. The intensities are also qualitatively correct with the lower energy-transfer 
peaks being less intense than the higher energy-transfer peaks. The gap is too small in our 
calculation, however, this could be accounted for via a GW gap correction. Note that we 
have not included the elastic scattering contribution in our calculation of the valence RIXS, 
which could effect the overall asymmetry of the signal. Finally, for all three spectra, the 
main edge occurs dbt cL Icir ger energy in the calculated results than in the experiment. This 
could be due to strong correlation effects, which would be expected to shift the Ti ci-states 
closer to the p-states. Another possible explanation which is important in the case of Ti K 
pre-edge XAS of Rutile Ti02 is the failure of the spherical muffin-tin approximation.— 
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B. Lifetime Broadening Suppressed XAS 



In addition to the RIXS planes, there are also several methods for obtaining life- 
time broadening suppressed (LBS) XAS. In high energy resolution fluorescence detected 
(HERFD) XAS,— an approximate absorption spectrum is found by partial fluorescence 
yield using a detector with resolution higher than the natural width due to core-hole life- 
time effects. This corresponds to viewing the spectra of constant emission energy in the 
RIXS plane. Another method of obtaining LBS XAS is to set the incident energy at a point 
well below the edge while scanning the emission energy.— Under certain assumptions, the 
spectrum obtained in this way is approximately proportional to the XAS signal multiplied 
by a Lorentzian with the width of the intermediate state core-hole. In Fig. IHlwe show a com- 
parison of experimental Cr-K edge HERFD XAS of K2Cr04 with our calculated results. We 
find reasonable qualitative agreement, with the exception of the peak just above 6000 eV, 
which is not seen in the calculation. We note however, that the size of this peak is sensitive 
to distortions. In addition, the amplitude of the main peak after the rising edge is too small. 
This could be due to the approximate treatment of the core-hole interaction or corrections 
to the spherical muffin-tin potentials used in FEFF9. Fig. H] shows a comparison of bulk 
metallic Ft L3 edge XANES compared to HERFD XAS and our calculated results. Here we 
see that all of the features are well reproduced although the broadening is too large at high 
energies, and the higher energy peaks are also red shifted toward the edge in comparison to 
the experimental result. 



V. CONCLUSIONS 



We have presented a theory of resonant inelastic x-ray scattering (RIXS) which is 
amenable to practical calculations as an extension of current x-ray-absorption and -emission 
codes. Starting from the Kramers-Heisenberg equation, we derive an expression for the 
RIXS cross-section which can be calculated using the real-space Green's function approach 
in the FEFF9 code. Inelastic losses and quasi-particle effects are included in terms of an 
effective spectral-function that is obtained from a quasi-boson model Hamiltonian. These 
many-body effects are incorporated into a single-particle approximation via a convolution 
with an effective spectral function. Quasi-particle self-energy effects are included based on 
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FIG. 3: (color online) Experimental (crosses) Cr K-edge HERFD XAS'*'^ of K2Cr04 compared to 
our calculated results (solid). 

a many-pole model of the dielectric function. Approximation of the many-body states as a 
product of an N-2 electron state with either two core electronic states (i.e. the ground state), 
or a core and photo-electron state (intermediate and final states) gives the cross section in 
terms of effective single particle Green's functions. The further approximation that the in- 
termediate and final photo-electron states are orthogonal with identical energies (valid at 
high photo-electron energies) gives the signal in terms of a convolution of the XAS and XES 
spectra. In addition, we have derived a formulation of the core-core RIXS spectrum that, due 
to the localized nature of AV, depends primarily on the Green's functions evaluated close 
to the absorbing atom. The extent of this localization is yet to be thoroughly investigated, 
although the degree of agreement between our calculations and experimental data suggests 
that the on-site approximation is valid for the systems shown here. In addition, the on-site 
approximation provides good qualitative agreement for the core- valence RIXS, although the 
approximation is less justifiable. The theory is implemented in an efficient program which is 
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FIG. 4: (color online) Pt L3 edge normal XANES (crosses) compared with HERFD XAS (x's) and 
om' calculated result (solid). 

an extension of the FEFF9 real-space muhiple-scattering code, which calculates the RIXS 
spectrum for core-core as well as core- valence RIXS. Several illustrative calculations were 
presented within the quasi-particle approximation where the spectral function is replaced 
by a 5-function, which appears to be a reasonable approximation for these cases. Calcu- 
lated results for MnO and for Anatase Ti02 based on this quasi-particle approximation are 
found to agree qualitatively with experimental spectra: the results reproduce both pre-edge 
and main edge features, the behavior of pre-edge features with varying core-hole interaction 
strength, and peak structure due to solid state effects in valence RIXS. Further investigations 
including treatments beyond the quasi-particle approximation of Eq. (!33|) will be reserved 
for the future. It should be noted that our current formalism for the spectral function is 



not suited for highly correlated materials, a. 
Hamiltonian is possible as suggested in Ref. 



though an extension of the quasi-boson model 
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Appendix A: Local behavior of RIXS: Derivation of T{D,) 

The RIXS cross section is given in terms of a product of three Green's functions, i.e., 
lm[g^{Ei)g'^{E2)g^{Eiy]. However, this expression is not very useful, since the spacial ar- 
guments of the Green's function must be integrated over all space to obtain the resonance 
l/(£'l — E2). In order to sec this, imagine that we ignore the core hole potentials in both 
Green's functions. In this case, we may write the Green's functions in spectal representa- 
tion, and noting that the wave functions are now orthonormal, we may rewrite the above 
expression as \lm[g{E2)/{El — E2 + iF)]^, which is what we expect for energies well above 
threshold, where the effect of the corehole potential is neghgable. While this hmit is easy 
to show analytically, doing so numerically within the real-space MS Greeen's function for- 
malism proves quite difficult. Below, we derive an alternative expression for the RIXS cross 
section which takes advantage of the localization of core-hole potential. Let us first define 
the Hamiltonian operators corresponding to the deep (b) and shallow (c) core holes. 

hc^ho + V, (Al) 

h = ho + Vb = K + AV (A2) 

Ay = H - K- (A3) 

We can use these definitions to rewrite the Green's funtions as follows 

//le=-l + £^2/, (A5) 

where we have left the energy arguments off of the Green's functions for the sake of brevity. 
Using these relations gives 

, , ^ g''-g' + g'AVg'^ ^ D^g' - g' 

9 9 TP T? T7< TP ^ ^ ) 

y y 771* 771 771* TP '' ^ ^ 
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where D = 1 + AVg'^^. Applying the above relations to g^{Ei)g'^{E2)g'^{Eiy gives 



Wglg'^ + g'ig^'g'^)] 

{D^g--g^)g^^ , g\g-D-h^^) 



El — E2 
D^g^D 1 



El - E2 



, gV^ , ^ 



\Ei — -E^2P El — i?; 



Noting that the second term above is real, we have 



(A8) 



(A9) 



Finally, we see that the transition matrix element T defined in Section (IIT]) can be related 

T = DPd (AlO) 



to D, i.e.. 



Appendix B: Effective Spectral function 



As shown in Sec. ([IT]), the RIXS cross section is given in terms of the ground state 
expectation value of a product of three many-body Green's functions. Here we will derive 
an expression based on quasi-particle Green's functions and a many-body spectral function. 

{^o\G\Ei)G''iE2)G''\Ei)\^o) 

nin2 

X {KjG'\Ei)e-''\^l) (Bl) 

Note that if we are expanding to second order in the boson couplings, only |$„) containing 
single boson excitations, = aJ,|$o) contribute. Now, 

1$^) = Z^e-^^l*^) (B2) 

where = [a^^ + AV/un] and AV" = V^"" - V^. Thus Eq. ^ becomes 

{%\G\Ei)G'^{E2)G'^\Ei)\%) 

nin2 

x«|G^ni5i)e-">^). (B3) 
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We now define the amplitudes to create a single boson as 

yn AV" 

a\ = V'd'iE^) < = V^G^E^) , (B4) 

and the amplitudes to annihilate a boson as 

yn AV"' 

/3^ = G''(Ei)(V^")*-^; /3^ = 6-^(^2) (y")* . (B5) 

Then combining Eq. (B3)-(B5) gives six terms quadratic in these amplitudes plus the zeroth 
order term. If wc neglect the ofT-rcsonant terms, i.e. those terms which contain a G^, G**^ 
with difi^erent energy arguments, this leaves only two terms plus the zeroth order term, 

{<^>o\G\E,)G''{E2)G'\E,)\%) 

^e-'^'i^G',{E^)G^,{E,)G'^iE^) 

+ G',(E,)P:Gl(E2)alG'^(E,) 

+ PlGl{E,)G'i,{E,) , (B6) 

where Gq — ($o|G*|$o), and G\ — This result can now be written in terms of 

a double convolution with a spectral function, i.e., 

{^o\G\E^)G%E2)G''\E^)\^^) 

= j dwidu-2 A^s{Ei,E2,u:i,uj2) 

X Gl{E^ - u{)Gl{E2 - uj2)G'J{E, - u,), (B7) 
where the effective spectral function is given by 

Aes{Ei, E2,UJi,UJ2) 

+ J2 mE2hn{E2)SMS{u2 - 

n 

+ \/3!^{E,)\'6{ui - UJn)S{uJ2 - UJn)] \ • (B8) 
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